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i ABSTRACT 

We present near-infrared spectra of the nuclear disk in the nearby radio 
galaxy NGC 5128 (Centaurus A). On the basis of the observed strengths of 
the [S III] A0.9532/im and [Fe II] A1.2567/im lines, we classify NGC 5128 as a 
LINER. Modeling of the strengths of these and additional lines suggests that 
the nuclear region is powered by shocks rather than photoionization. 
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The galaxy NGC 5128 hosts the extragalactic radio source Centaurus A. It is our 
nearest radio galaxy, and has been the subject of intense study which has revealed a jet 
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Introduction 



seen at radio, near-infrared, and X-ray wavelengths (Meier et al. 1989; Joy et al. 1991; 
Schreier et al. 1979), and an unresolved nucleus which appears longward of 2 /im (Turner et 
al. 1992). Perhaps the most obvious feature of NGC 5128, however, is the dust lane which 
straddles the galaxy. Such prominent and well-ordered dust obscuration is more commonly 
associated with spirals, and is unusual for an early-type galaxy (although perhaps not for a 
radio galaxy: de Koff et al. 1996). It is widely believed that the dust lane was formed from 
material stripped from a dusty disk galaxy (Baade & Minkowski 1954; Thomson 1992), and 
a natural consequence of this merger would be substantial star formation, which is indeed 
inferred (Quillen, Graham & Frogel 1993; Storchi-Bergmann et al. 1997). 

The type of spectroscopic study normally performed on the nuclei of active galaxies, 
however, has yet to be undertaken on NGC 5128. Since the nucleus is not visible in 
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ground-based optical images, no astrometry has been performed to allow it to be acquired 
for spectroscopy. Phillips (1982) assumed the nucleus was located at the near-infrared 'hot 
spot' of Kunkel & Bradt (1971), a feature which can now be identified with the infrared 
jet, and Storchi-Bergmann et al. (1997) used a wide (10") slit which contained a significant 
contribution from material excited by the young stars in the larger-scale disk. Recently, 
Schreier et al. (1998) have shown that there is a physically distinct Paa-emitting region 
at the nucleus approximately 2" in size, which they interpret as a disk. As similar, albeit 
larger, disks have been seen in a number of FR I radio galaxies, a nuclear spectrum might 
reveal whether this smaller-scale disk shares their properties. 

The only spectrum of the nucleus alone is that of Meadows & Allen (1992, MA92), who 
used the Infrared Imaging Spectrometer on the 3.9-m Anglo-Australian Telescope, IRIS 
(Allen et al. 1993), which has the capability to image a field before inserting the dispersing 
element. In this paper, we present further spectra taken with IRIS, covering the wavelength 
range 0.9 /iin < A < 2.5 /zm, and use them to investigate the excitation mechanism operating 
at the nucleus of NGC 5128. 

Throughout this paper, we adopt a distance of 3.6 Mpc to NGC 5128 (Soria et al. 
1996), corresponding to a projected linear scale of 17pc arcsec -1 . 

2. Observations and reduction 

Data were taken with both narrow (I'M) and wide (5"8) slits at a spatial resolution of 
0"8 pixel -1 . IRIS was configured to use the IJ (0.9-1.5 /zm) and HK (1.4-2.5 /jm) echelles 
and provided a spectral resolution of A/AA ~ 400 (narrow slit) or ~ 100 (wide slit). The 
true nucleus of NGC 5128 was acquired by locating the continuum peak at K', and spectra 
were taken at either end of the 13" slit and differenced to remove the sky and underlying 
galaxy. All data were taken with the slit oriented east-west in approximately 2" seeing. 
Narrow slit HK echelle data were taken on the nights of UT 1991 April 26/27 and 1991 
June 25/26 for a total integration time of 7400 s, and were presented in MA92. Narrow slit 
IJ data totaling 2000 s were taken on UT 1992 April 14. Wide slit data were taken on UT 
1992 September 13, with an exposure time of 1200 s in each echelle. 

Standard spectral reduction techniques were applied to remove the spatial and spectral 
response of the detector, the curvature of the echelle orders, and the apparent rotation of 
the slit, which changes across the orders, before one-dimensional spectra were extracted. 
The narrow slit spectra, which were taken in non-photometric conditions, have been flux 
calibrated by scaling the continuum fluxes so that they matched photometry from the 



- 3- 



images of MA92, which in turn were calibrated using the aperture photometry of Turner et 
al. (1992). We estimate that the systematic uncertainties introduced by this method are 
less than 20%; the ratios derived from line pairs within the IJ or HK echelles are, of course, 
unaffected. Observations of standard stars were used to calibrate the wide slit data. 



3. Results and analysis 

The narrow slit spectrum is shown in Figure |l|. We present the extracted line fluxes 
from both the wide and narrow slit spectra in Table |TJ. These fluxes were obtained by 
measuring the signal above a linear continuum level which was determined by eye. The 
placement of this continuum is the dominant source of uncertainty in the measurements, 
and so several measurements with different 'reasonable' continuum levels were made to 
estimate the possible errors. The measurement uncertainties from the narrow slit data are 
much smaller than those from the wide slit spectra, due to the lower sky noise, which allows 
more accurate continuum determination. 

Given the poor seeing conditions of our observations, we would expect the line fluxes 
from the wide slit spectra to be ~ 40% higher than those from the narrow slit spectra due 
to slit losses, even if the line emission were unresolved. Coupled with the uncertainty in the 
flux calibration of the narrow slit data, we are unable to find any evidence for significant 
spatial extension in any of the emission lines. It therefore seems likely that all the emission 
lines share the morphology of Paa seen by Schreier et al. (1998). 

3.1. Reddening to the emission line region 

The most reliable estimate of the reddening to the emission line region comes from the 
near-infrared [Fe II] lines, whose intrinsic intensity ratio is fixed (Nussbaumber Sz Storey 
1988). The observed ratio in the narrow slit data suggest Ay = 3 mag, although it is also 
consistent with zero reddening, and provides a formal 3a upper limit of Ay < 8. In the 
wide slit data, the 1.6435 //m line is strongly affected by poor sky subtraction and so no 
useful ratio can be obtained. 

The only reliable H I line ratio we have from the wide slit data is Pa/3/Br7 = 1.25 ±0.25. 
Compared to the Case B ratio of 5.88, this implies an extinction of Ay = 10.9 ± 2.2. 
Although the Pa/3/Br7 and Pa7/Br7 ratios from the narrow slit data are less reliable due 
to the method of flux calibration, they too support a value for the extinction of about ten 
magnitudes. However, the observed narrow slit Pa/3/Pa7 = 1.83 ± 0.35 is consistent with 
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the Case B value of 1.81, implying negligible reddening, whereas if Ay ph 10, it should be 
twice as large. A formal 3a upper limit to the extinction based on this ratio (which is not 
affected by flux calibration uncertainties) is Ay < 7. 7 mag. We therefore conclude that the 
Br7 flux is enhanced by an additional component which, presumably by virtue of being 
heavily reddened, does not contribute to the Paschen lines. This component could arise 
from a putative broad line region. Assuming the broad line region suffers the Ay ~ 30 mag 
of extinction which obscures the nuclear continuum (Giles 1986; MA92) then any broad 
Br7 will be reduced in flux by a factor of about 20. If the broad line is intrinsically 40 
times brighter than the narrow line (e.g. Jackson & Eracleous 1995), then the ratio of 
narrow Pa/3/Br7 rj 4, corresponding to Ay ph 3, in line with that derived from the [Fe II] 
ratio. This extinction is also consistent with obscuration from the dust lane alone, based 
on the observed nuclear J — H color if we assume that approximately half the nuclear H 
band emission is non-stellar (MA92), and so we adopt Ay = 3 mag for the extinction to 
the emission line region. This value is somewhat lower than that derived by Schreier et 
al. (1996), but by failing to correct for the presence of the heavily-reddened nucleus these 
authors are likely to have overestimated the extinction to the NLR. 



3.2. Excitation mechanism 

MA92 suggested a starburst origin for the emission lines on the basis of their HK 
spectrum, and in particular the H2 v =1-0 S(l)/Bry ratio. However, the low dissociation 
temperature of H 2 means the H 2 line emission is not produced cospatially with the other 
emission lines, and this line ratio is therefore not ideal for determining the excitation 
mechanism. The observed ratio of 0.6 is equally consistent with starburst, Seyfert, and 
LINER galaxies (Moorwood Sz Oliva 1988; Larkin et al. 1998). A starburst can, however, be 
excluded on the basis of the luminosity of the [Fe II] emission. Following the discussion of 
van der Werf et al. (1993), the radius of an individual SNR at the beginning of the radiative 
phase is ~ 6 pc, and our spectroscopic observations indicate that the [Fe II] emission is at 
best barely resolved. Even making the generous assumption that the line-emitting region 
is spherical with a diameter of 20 pc, there is only room for at most 5 SNRs but, adopting 
L[Feii] = 10 3 1/© for a single remnant, about 20 are needed to provide the extinction-corrected 
luminosity. 

Of the emission lines we observe, those of [Fe II] and [S III] have been the most 
well-studied, and we present a line ratio diagram in Figure |2|. Although few objects have 
been observed in both emission lines, resulting in a sparsely populated plot, many more have 
fluxes for one or other of the two lines, and different classes of object are fairly well-separated 



5 



in each of the ratios, although there are regions of overlap. We can therefore subdivide the 
parameter space as illustrated. We also plot the theoretical line ratios as determined using 
the photoionization code cloudy V90.04 (Ferland 1996) for a plane-parallel slab with 
densities in the range tt-h = 10 3 -10 6 cm -3 illuminated by a power law spectrum extending 
from 10 /im to 50keV with S v oc u~ 1A , a range of ionization parameter U = 10~ 2 -10~ 5 and 
a gas phase iron abundance one-twentieth of solar. The observed ratios can be explained by 
an ionization parameter U ~ 10 -4 , typical of LINERs, (e.g. Ferland & Netzer 1983), and 
a density hh ~ 10 4 ' 5 cm -3 . This relatively high density is required to produce the H I line 
emission, since the extinction-corrected Ha luminosity is L-^ a = (2.2 ±0.3) x 10 32 W (under 
Case B), which is larger than the luminosity in the nuclear disk of M 87 (Ford et al. 1994), 
yet the emission in NGC 5128 comes from a much smaller region. 

We observe He I A1.0830/xm/Pa7 ~ 4. This compares with typical values of ~ 2 (from 
Case B recombination and the observed He I AO. 4471 /xm/H/3 ratio) in both H II regions 
(Doyon, Puxley & Joseph 1992) and AGN (Ferland & Osterbrock 1986), although larger 
ratios of 6.2 and 6.7 have been observed in Orion and NGC 4151, respectively (Osterbrock, 
Shaw & Veilleux 1990), and a ratio of 5.8 is predicted for the starburst galaxy NGC 3256 
based on the observed He I A2.06 yum/Br7 ratio (Doyon, Joseph & Wright 1994). This ratio 
is therefore unable to tell us anything about the excitation mechanism. 

We have also detected the [C I] AA0.9840,0.9850 /im and [S II] A1.0330/xm blends. They 
have also been seen in NGC 4151, where their fluxes relative to [S III] AO. 9532 /im were 0.01 
and 0.05, respectively (Osterbrock et al. 1990). We would expect stronger emission from 
these low-ionization lines in a source with a lower ionization parameter, but there are no 
data from such sources with which to make comparisons, so we must compare our data with 
model predictions. We again use cloudy for the photoionization models and mappings II 
(Sutherland & Dopita 1993) to investigate the effects of shock excitation, performing the 
same self-consistent iterative process as Dopita & Sutherland (1996). We used a solar iron 
abundance for the shock models, since grains would be destroyed by the shocks, liberating 
iron into the gas phase (Greenhouse et al. 1991). A high preshock density is required, since 
the scaling relation for H/3 luminosity of Dopita & Sutherland (1996) for a planar shock 
with speed V covering an area d 2 , requires 



to explain the inferred H/3 luminosity. 

Table |] presents the dereddened line fluxes and model results (relative to Pa/3). As 
mappings II does not provide the strength of the He I A1.0830 /im line in its output, we are 
unable to compare the observed value with its prediction. Instead we present a lower limit 
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to the strength of this line, from Case B ratios and the strength of He I AO. 4471 /zm. The 
relatively high density required will result in a somewhat larger flux, since He I A1.0830/im 
is readily enhanced by collisional excitation. Since all the emission lines we are studying 
arise from a fairly small wavelength region, the line ratios are not strongly affected by 
reddening. Even if the true extinction were as high as Ay = 5 mag, the most strongly 
affected line ratio, [S III] /Pa/3, would only increase by 40%, not altering the conclusions we 
draw. 

Although we have not investigated all the parameter space for the shock models, due 
to the computational effort involved, it is clear that shocks provide a rather better fit to the 
observed line ratios than do the photoionization models. In particular, photoionization fails 
to produce the observed line ratios over the range of ionization observed (C°, S + , S 2+ ), as 
has often been seen in the ultraviolet (e.g. Dopita et al. 1997). We therefore favor shocks 
over photoionization, although ideally we would like some empirical evidence to support 
this conclusion, from observations of the same emission lines in objects where shocks are 
known to be the dominant excitation mechanism. 



4. Discussion 

The above arguments support the classification of NGC 5128 as a LINER, even though 
the emission lines which define this class (Heckman 1980) are not observable. While 
Storchi-Bergmann et al. (1997) also made such a classification, their much larger aperture 
was dominated by emission lines originating in the starburst disk, and so did not directly 
inform on the nuclear properties. The fact that we see a shock-excited LINER spectrum in 
the nuclear disk is perhaps not surprising, since the same is seen in M 87 (Virgo A; Harms 
et al. 1994; Dopita et al. 1997) and the similarities in the radio properties of NGC 5128 and 
M 87 are quite extensive (Tingay et al. 1998). Whether shock-excited nuclear disks are a 
common feature of FR I radio galaxies is a question that requires further study. 

Spectroscopic observations in the ultraviolet might also confirm our claim that 
the emission lines are powered by shocks, since the predictions of shock excitation and 
photoionization for lines such as C IV AO. 1550 /im and [Ne V] A0.3426yum can differ by 
factors of 10 4 (e.g. Dopita et al. 1997). Although there is significant extinction to the 
emission-line region, the emission lines are sufficiently bright that these UV lines might 
still be detectable. Indeed, our shock model, coupled with Ay = 3 mag of extinction, 
predicts /civ ~ 2 x 10~ 18 Wm~ 2 , which is observable with the Space Telescope Imaging 
Spectrograph (STIS). However, we note that the C IV flux will be an order of magnitude 
fainter for every magnitude of visual extinction by which we have underestimated the 
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true obscuration, so while its detection would unequivocally rule out photoionization, its 
non-detection would not provide a definitive result. 

5. Summary 

We have presented near-infrared spectra of the nucleus of NGC 5128, covering the 
wavelength range 0.9-2.5 /im. The emission line ratios are characteristic of a LINER, 
and can be modeled fairly well by shock excitation or photoionization by a power law 
continuum although shocks provide a noticeably better fit. Coupled with the similarities 
between Cen A and Vir A, we favor shocks as the dominant excitation mechanism in 
NGC 5128. Ultraviolet spectroscopy might be able to confirm this through the detection 
of high-ionization lines (e.g. C IV AO. 1550 /im), but the high obscuration in the UV could 
easily push them below a feasible detection threshold. 

This work was performed at the Jet Propulsion Laboratory, California Institute of 
Technology, under a contract with the National Aeronautics and Space Adminstration. VM 
was supported for these observations by an Australian Postgraduate Research Award. We 
thank the Anglo- Australian Observatory for generous access to its computing and other 
facilities. We also wish to thank David Allen for assistance in acquiring the spectra and for 
many useful discussions, and the anonymous referee for constructive suggestions. 
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Table 1: Emission line fluxes through the different spectroscopic apertures. The I'M x 1'.'6 
aperture measurement have been calibrated as described in the text, and are subject to a 
~ 20% systematic uncertainty, in addition to the random errors listed. Fluxes indicated 
with a dagger (j) are strongly affected by an uncertain background level as a result of poor 
OH line subtraction. 
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Table 2: Comparison of excitation models. All fluxes have been normalized to Pa/5 = 100, 
and the observations have been corrected for reddening. The photoionization model has 
n H = 10 5 cm~ 3 and U = 10~ 4 . The shock excitation model is for a 120 km s" 1 shock with 
precursor in a medium with density nu = 3 x 10 4 cm" 3 and an equipartition magnetic field. 
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Fig. 1. — Nuclear spectrum of NGC 5128, extracted through a I'M x 1'.'6 aperture. The 
emission lines are labeled, as are regions of CO and atmospheric absorption. This spectrum 
was taken in non-photometric conditions, and the flux scale has been bootstrapped as 
described in the text. 
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Fig. 2. — [Fe II] /Pa/5 vs [S III] /Pa/9 line ratio diagram. The regions expected to be occupied 
by different classes of object are separated by dotted lines as shown in the inset, and those few 
objects for which data exist for both forbidden lines are plotted. We have plotted reddening- 
corrected line ratios for Seyfert galaxies (solid squares), starburst galaxies (open circles), 
and LINERs (open squares), using data from Kirhakos & Phillips (1989), Simpson et al. 
(1996), and Alonso-Herrero et al. (1997) and references therein, assuming Case B ratios to 
convert Ha and Br7 fluxes to Pa/5. The six-pointed star is the Galactic supernova remnant 
RCW 103 (Dennefeld 1986). The solid lines are theoretical predictions of photoionization 
models of different densities, as described in the text. 



